Exploitation of symmetry in periodic Self-Consistent-Field ab initio calculations: application to large three-dimensional compounds

Symmetry can dramatically reduce the computational cost (running time and memory allocation) of Self-Consistent-Field ab initio calculations for crystalline systems. Crucial for running time is use of symmetry in the evaluation of one- and two-electron integrals, diagonalization of the Fock matrix at selected points in reciprocal space, reconstruction of the density matrix. As regards memory allocation, full square matrices (overlap, Fock and density) in the Atomic Orbital (AO) basis are avoided and a direct transformation from the packed AO to the SACO (Symmetry Adapted Crystalline Orbital) basis is performed, so that the largest matrix to be handled has the size of the largest sub-block in the latter basis. We here illustrate the effectiveness of this scheme, following recent advancements in the CRYSTAL code, concerning memory allocation and direct basis set transformation. Quantitative examples are given for large unit cell systems, such as zeolites (all-silica faujasite and silicalite MFI) and garnets (pyrope). It is shown that the full SCF of 3D systems containing up to 576 atoms and 11136 Atomic Orbitals in the cell can be run with a hybrid functional on a single core PC with 500 MB RAM in about 8 h. © 2014 Science China Press and Springer-Verlag Berlin Heidelberg

Communication: Hole localization in Al-doped quartz SiO2 within ab initio hybrid-functional DFT

We investigate the long-standing problem of the hole localization at the Al impurity in quartz SiO$_2$, using a relatively recent DFT hybrid-functional method in which the exchange fraction is obtained \emph{ab initio}, based on an analogy with the static many-body COHSEX approximation to the electron self-energy. As the amount of the admixed exact exchange in hybrid functionals has been shown to be determinant for properly capturing the hole localization, this problem constitutes a prototypical benchmark for the accuracy of the method, allowing one to assess to what extent self-interaction effects are avoided. We obtain good results in terms of description of the charge localization and structural distortion around the Al center, improving with respect to the more popular B3LYP hybrid-functional approach. We also discuss the accuracy of computed hyperfine parameters, by comparison with previous calculations based on other self-interaction-free methods, as well as experimental values. We discuss and rationalize the limitations of our approach in computing defect-related excitation energies in low-dielectric-constant insulators.Comment: Accepted for publication in J. Chem. Phys. (Communications

5D quantum dynamics of the H2@SWNT system: Quantitative study of the rotational-translational coupling

The dynamics of the dihydrogen molecule when confined in carbon nanotubes with different chiral- ities and diameters are studied by using a 5 dimensional model considering the most relevant degrees of freedom of the system. The nuclear eigenstates are calculated for an (8,0) and a (5,0) carbon nanotubes by the State-Average Multiconfigurational Time-dependent Hartree, and then studied using qualitative tools (mapping of the total wave functions onto given subspaces) and more rigorous analysis (different kinds of overlaps with reference functions). The qualitative analysis is seen to fail due to a strong coupling between the internal and translational degrees of freedom. Using more accurate tools allows us to gain a deeper insight into the behaviour of confined species

Assessing Thermochemical Properties of Materials through Ab initio Quantum-mechanical Methods: The Case of α-Al2O3

The thermochemical behavior of α-Al2O3 corundum in the whole temperature range 0–2317 K (melting point) and under pressures up to 12 GPa is predicted by applying ab initio methods based on the density functional theory (DFT), the use of a local basis set and periodic-boundary conditions. Thermodynamic properties are treated both within and beyond the harmonic approximation to the lattice potential. In particular, a recent implementation of the quasi-harmonic approximation, in the Crystal program, is here shown to provide a reliable description of the thermal expansion coefficient, entropy, constant-volume and constant-pressure specific heats, and temperature dependence of the bulk modulus, nearly up to the corundum melting temperature. This is a remarkable outcome suggesting α-Al2O3 to be an almost perfect quasi-harmonic crystal. The effect of using different computational parameters and DFT functionals belonging to different levels of approximations on the accuracy of the thermal properties is tested, providing a reference for further studies involving alumina polymorphs and, more generally, quasi-ionic minerals

A new structural model for disorder in vaterite from first-principles calculations

Both of the previously proposed Pbnm and P6522 ordered structures for vaterite are found to be unstable transition states using first principles methods. Five stable structures are located, the lowest energy one being of P3221 symmetry. Since interconversion between these structures requires only thermal energy, this provides an additional source of disorder within the vaterite structure

Probing vibrational modes in silica glass using inelastic neutron scattering with mass contrast

The effective vibrational density of states (VDOS) has been derived from inelastic neutron-scattering data, for isotopically substituted Si O 18 2 and Si O 16 2 glasses, to gain information about the relative contribution to the Si and O partial VDOS. This is a necessary point of comparison for vibrational mode analyses of molecular-dynamics models. The mass contrast has led to a measurable shift between vibrational mode frequencies in the effective VDOS of Si O 18 2 and Si O 16 2, which is well reproduced in an ab initio simulation. The vibrational band centered at 100.2 meV is confirmed to have significantly lower contribution from the oxygen partial VDOS, than the higher (150.3 and 135.8 meV) and lower energy bands (53.3 meV)

Bulk and surface energetics of lithium hydride crystal: benchmarks from quantum Monte Carlo and quantum chemistry

We show how accurate benchmark values of the surface formation energy of crystalline lithium hydride can be computed by the complementary techniques of quantum Monte Carlo (QMC) and wavefunction-based molecular quantum chemistry. To demonstrate the high accuracy of the QMC techniques, we present a detailed study of the energetics of the bulk LiH crystal, using both pseudopotential and all-electron approaches. We show that the equilibrium lattice parameter agrees with experiment to within 0.03 %, which is around the experimental uncertainty, and the cohesive energy agrees to within around 10 meV per formula unit. QMC in periodic slab geometry is used to compute the formation energy of the LiH (001) surface, and we show that the value can be accurately converged with respect to slab thickness and other technical parameters. The quantum chemistry calculations build on the recently developed hierarchical scheme for computing the correlation energy of a crystal to high precision. We show that the hierarchical scheme allows the accurate calculation of the surface formation energy, and we present results that are well converged with respect to basis set and with respect to the level of correlation treatment. The QMC and hierarchical results for the surface formation energy agree to within about 1 %.Comment: 16 pages, 4 figure

Elastic and magnetic effects on the infrared phonon spectra of MnF2

We measured the temperature dependent infrared reflectivity spectra of MnF2 between 4 K and room temperature. We show that the phonon spectrum undergoes a strong renormalization at TN. The ab-initio calculation we performed on this compound accurately predict the magnitude and the direction of the phonon parameters changes across the antiferromagnetic transition, showing that they are mainly induced by the magnetic order. In this material, we found that the dielectric constant is mostly from phonon origin. The large change in the lattice parameters with temperature seen by X-ray diffraction as well as the A2u phonon softening below TN indicate that magnetic order induced distortions in MnF2 are compatible with the ferroelectric instabilities observed in TiO2, FeF2 and other rutile-type fluorides. This study also shows the anomalous temperature evolution of the lower energy Eu mode in the paramagnetic phase, which can be compared to that of the B1g one seen by Raman spectroscopy in many isostructural materials. This was interpreted as being a precursor of a phase transition from rutile to CaCl2 structure which was observed under pressure in ZnF2.Comment: 8 pages, 8 figures, updated version accepted in PR